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ABSTRACT 

The conversion of theoretical, as well as geometric heights to practical heights requires the application of 
geoidal undulations from a geoid model. The various global geopotential models that are readily available for 
application in any part of the world do not best-fit regions, as well as countries. As a result, there is a need to 
determine the local geoid models of local areas, regions and countries. This study determines the local geoid 
model of Kampala in Uganda for orthometric heights computation by comparing three plane geometric geoid 
surfaces. A total of 19 points were used in the study. The least squares adjustment technique was applied to 
compute the models’ parameters. Microsoft Excel programs were developed for the application of the models 
in the study area. The Root Mean Square Index was applied to compute the accuracy of the models. The three 
geometric geoid models were compared using their accuracy to determine which of them is most suitable for 
application in the study area. The comparison results show that the three models can be applied in the study 
area with more reliability, with greater confidence in model 2. 
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INTRODUCTION 

Research and Development of a system of Global Navigation 
Satellite System (GNSS) positioning have made a method of 
surveying faster, easier and more accurate than before in the 
application of both scientific and applied geodesy and 
especially for geoid modelling when gravity data is 
unavailable. The geoid according to Heiskanen and Moritz 
(1967) is the “mathematical figure of the earth”. It is the 
surface which coincides with the mean sea level assuming that 
the sea was free to flow under the land in small frictionless 
channels. Ubajekwe (2011) explained further that the mean sea 
level is not quite an equipotential surface owing to non- 
gravitational forces (such as ocean currents, winds and 
barometric pressure variation). In geodetic surveying, the 
computation of the geodetic coordinates of points is commonly 
performed on a reference ellipsoid closely approximating the 
size and shape of the earth in the area of the survey. The actual 
measurements made on the surface of the earth with certain 
instruments are however referred to the geoid (Oduyebo et al., 
2019). The geoid may be obtained by modelling from among 
others to serve as a vertical reference, as well as a datum for 
height determination. A global geoid is designed as global best- 
fit for the whole world with each country adopting her own 
local or regional version to fit her own mapping needs and 
purposes, as well as serve as a datum in engineering 
specifications. 

The geoid surface being a continuous equipotential surface is 
used in defining the heights vertical reference in place of 
unreliable mean sea level. Local geoid modelling of 
areas/regions has become crucial as the GNSS ellipsoidal 
heights are not practical heights used in engineering 


constructions. The GNSS ellipsoidal heights are theoretical 
heights obtained by computation on a specified ellipsoid. To 
apply the said heights in engineering constructions, a 
conversion is necessary. To do this, a geoid model is required 
to obtain the occupied points/controls respective geoid heights. 
Having obtained the geoid heights of the points from the geoid 
model, the orthometric heights of the respective points are 
obtained using the relation (Oluyori, at al. 2018): 

H = h-N 

(i) 

Where, 

H = Orthometric height 
h = Ellipsoidal height 
N = Geoid height 

Different studies have determined the local geoid models of 
various areas with different accuracy achieved. Oduyebo et al. 
(2019) determined the local geoid model of Benin City by 
comparing three gravimetric-geometric geoid models of the 
study area with an accuracy of 0.6746m. Oluyori et al. (2018) 
also determined the local geoid model of the Federal Capital 
Territory, Abuja, where the accuracy of 0.419m, was obtained 
by comparing GNSS/Levelling and EGM 2008 geoidal 
undulations. 

Although, Uganda gravimetric geoid model of 2014 
(UGG2014) had been experimented successfully by the KTH 
method for orthometric height determination using the GNSS 
equipment (Ssengendo et al., 2015a). The experiment was 
carried out using a global geopotential model (EIGEN-6C4). 
Global geopotential models are determined for application in 
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any part of the world. They best fit the entire globe but not 
local areas, regions or countries. Their application at local 
areas, regions and countries, yield less accurate results. 
Consequently, there is a need to determine local geoid models 
of areas, regions and countries. This study comparatively 
analyses three geometric geoid surfaces using observed GNSS 
ellipsoidal heights and existing orthometric heights to 
determine the most suitable geoid surface for application in the 
Kampala area. 

The geometric method of local geoid modelling is applied to 
small areas. A study has stated that a small area is the one that 
is less than 200km 2 in size (Schofield, 2007). In surface fitting 
methods, models adopted depend on the size and nature of the 
variation of the heights of the points used as opined Oluyori 
(2019). As such, Romans (2007) suggested that for small areas, 
models with four parameters or lower, be adopted for better 
accuracy. Kampala has an area of 196km 2 which by Romans 
(2007) was classified as a small area in this study. 

Objectives of the study are: to determine ellipsoidal heights of 
controls from dual-frequency DGPS observations; to obtain the 
geoidal heights (N) of the controls by finding the differences 
between the GNSS ellipsoidal heights and the existing 
orthometric heights; to develop Microsoft Excel programs to 
compute the models' parameters, geoid heights, orthometric 


heights and the accuracy, and to compare the computed 
models’ accuracy to obtain the most suitable model for 
application in the study area. 

The Study Area 

Kampala is the capital city of Uganda and occupies a series of 
hills at an average elevation of 1,190m and is located north of 
Lake Victoria. It has an area of 196km 2 . The City today has 
grown into a Commercial, Educational, Cultural and 
Administrative Centre of Uganda with an approximate 
population of 2.5 million people. Considering that the 
population was 330,700 in 1969, 1,208,544 in 2002, and 
1,811,794 in 2010, this signals rapid urbanization in the 
country (Oonyu and Esaete 2012). Quite a lot of survey and 
mapping, as well as engineering activities, must be expected. 
The study area falls within the UTM zone 36N rectangular 
coordinate system. 

Topography 

Kampala is made up of Central, Kawempe, Makindye, Nakawa 
and Lubaga divisions and it sits on hills including; Kasubi Hill, 
Mengo Hill, Kibuli Hill, Namirembe Hill, Rubaga Hill, 
Nsambya and old Kampala Hill. Figure 1 shows the location of 
Kampala in Uganda and its five divisions as located within its 
boundaries. 



Figure 1: Administrative Units of Kampala and its Location in Uganda 
Source: Irumba (2015) 
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Plane Geometric Geoid Surfaces 

Plane geometric geoid surfaces are geoid models with the highest degree as 1. They are used to model the local geoid of small 
areas. The plane geometric geoid surfaces given by Abdallah (2010) as model 1, model 2 and model 3 are respectively 

N t =a 0 -a l x i +a 2 y i (2) 

Y =a 0 + a x x t + a 2 y. + a 3 x t y. (3) 

N { = a 0 + a l x i + a 2 y t + a 3 Ah (4) 

Where, 

N f = Geoid height 

X- = Northing of the observed point 

y t = Easting of the observed point 

Ah = Difference between the average ellipsoidal heights of the observed points and individual point 
CLq , a 1 , d 2 , CL 3 = Model parameters 

The application of equations (2), (3) and (4) in local geometric geoid modelling, requires the use of least squares technique to 
obtain the models’ parameters. And it involves the writing of observation equations. Here, the number of observation equations to 
be written must equal the number of the observed points. The least squares models and the procedures used for the computation 
of the geometric geoid model parameters are detailed in Eteje and Oduyebo (2018) and Eteje, et al, (2018). 

Root Mean Square Error Computation (RMSE) 

The Root Mean Square Error (RMSE) is an indicator of accuracy. It is used for the computation of the accuracy of local 
geometric geoid models. Its application for accuracy computation in geoid modelling involves the comparison of the computed 
geoid heights obtained from the differences between the ellipsoidal and the orthometric heights and the model geoid heights of 
points. Also, the RMSE of the geometric geoid model can be obtained using the known and the model orthometric heights. The 
RMSE index used for accuracy computation as given by Eteje and Oduyebo (2018) is 


Where, 


RMSE = ± 




vfv 

n 


V = (N/H) known -(N/H) uodel 

(A/7 H) known = Point known geoid/orthometric height 

(N / H) MODEL — Point model geoid/orthometric height 
n = Number of points 


(5) 


The Gravity Field Model EIGEN-6C4 

According to Kostelecky et al (2015), EIGEN-6C4 (European 
Improved Gravity model of the Earth by New techniques) is a 
static global combined gravity field model up to degree and 
order 2190. It has been elaborated jointly by GFZ Potsdam and 
GRGS Toulouse and contains the following satellite and 
ground data: 

1. LAGEOS-1/2 (deg. 2 - 30): Satellite Laser Ranging data 
1985 - 2010 

2. GRACE, GNSS-SST and K-band range-rate data, 
processing according to RL03 GRGS (deg. 2 - 130): ten 
years 2003 - 2012 

3. GOCE, Satellite Gravity Gradiometry (SGG) data, 
processed by the direct approach including the gravity 
gradient components Txx, Tyy, Tzz and Txz out of the 
following time spans 837 days out of the nominal mission 
period 20091101 - 20120801, 422 days out of the lower 
orbit phase between 20120801 - 20131020. 

4. Terrestrial data (max degree 370): DTU12 ocean geoid data 
and an EGM2008 geoid height grid for the continents. The 


combination of these different satellites and surface data 
sets has been done by a band-limited combination of 
normal equations (to maximum degree/order 370), which 
are generated from observation equations for the spherical 
harmonic coefficients (Shako et al, 2013, Kostelecky et 
al, 2015). The resulted solution to degree/order 370 has 
been extended to degree/order 2190 by a block diagonal 
solution using the DTU10 global gravity anomaly data 
grid. 

The EIGEN6C4 has been evaluated and adopted for geodetic 
computation, as well as geoidal heights and orthometric heights 
determination in Uganda. 

METHODOLOGY 
Data Acquisition 

A total of 19 points, were used in the study, 12 points for the 
modelling of the geometric geoid of the study area, 4 points 
(U2004 to U2007) for validation using the existing global 
geoid model (EIGEN6C4) (MLHUD, 2019) (See Figure 2) and 
3 points (71Y121, 71Y126 and 71Y147) also for validation 
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using orthometric heights from spirit levelling. Points U2004 to 
U2007 are respectively located at Road reserve, Kawali Lweza; 
UNRA Road reserve, Kinaawa; Island of Kira road to 
Bukoto/Kamwokya, Bukoto and Makerere university business 
school, Nakawa. The points were observed using Trimble R7 
receivers to obtain their coordinates and ellipsoidal heights. 4 
of 19 of the existing/known orthometric heights of the points 
were obtained from the Ministry of Land, Housing and Urban 


Development (MLHUD) and 15 of 19 were obtained by a 
private firm SIG in 1993. The orthometric heights were also 
obtained from the MLHUD. They were computed using the 
ellipsoidal heights from GNSS observation carried out with 
Leica 1200 GPS receivers and accessories and the geoid 
heights obtained from a global geoid model (EIGEN6C4) 
adopted for geodetic computation in Uganda. 



Figure 2: Gravimetric Orthometric Heights Test Points 
Source: Ministry of Land, Housing and Urban Development (MLHUD) (2019) 


Data Processing 

The GNSS observations were processed in the UTM zone 36N 
on the GRS80 ellipsoid using Bernese (Version 5.0) to obtain 
the northing, easting and ellipsoidal heights of the occupied 
controls. The geoid heights of the 12 points used for the local 
geometric geoid modelling, were computed by finding the 
differences between the ellipsoidal heights from GNSS 
observation and the known orthometric heights of the points. 
Also, the geoid heights of the 3 points used for validation, were 
computed using the ellipsoidal heights from the GNSS 


observation and the orthometric heights from the spirit 
levelling. The geoid heights of the points were computed using 
equation (1). Those of the 4 validation points, were obtained 
from EIGEN6C4 using the coordinates of the points. The 
existing/known orthometric heights were also computed using 
equation (1). The levelling data of the 3 validation points, were 
reduced using the height of the instrument/collimation method. 
Table 1 shows the coordinates, ellipsoidal heights and the 
known orthometric heights of the 12 points used for the 
geometric geoid modelling. 


Table 1: Coordinates, Ellipsoidal Heights, Geoid Heights and Known Orthometric Heights of the 12 Points 

Geometric Geoid Modelling Points 


STATION 

EASTING 
(x) (m) 

NORTHING 

(y) (m) 

ELLIPSOIDAL 
HEIGHT, h (m) 

Known 

ORTHOMETRIC 
HEIGHT, H (m) 

GEOID HEIGHT, 
N (m) 

71Y65 

449421.520 

39563.190 

1209.921 

1222.122 

-12.215 

71Y80 

451601.150 

37537.100 

1253.079 

1265.410 

-12.334 

71Y97 

451574.520 

37388.890 

1255.641 

1267.942 

-12.337 

71Y125 

454888.680 

34774.920 

1137.505 

1150.081 

-12.460 

71Y141 

457725.480 

30265.860 

1176.309 

1188.934 

-12.576 

71Y143 

457490.680 

26773.960 

1152.722 

1165.290 

-12.639 

71Y149 

458258.050 

35349.840 

1151.609 

1164.099 

-12.531 

71Y151 

458861.180 

39178.510 

1237.155 

1249.643 

-12.491 

71Y152 

458281.140 

40110.200 

1211.352 

1223.836 

-12.443 

71Y153 

458393.460 

42176.260 

1193.261 

1205.659 

-12.403 

71Y154 

457752.180 

44442.970 

1174.250 

1186.604 

-12.365 

71Y155 

459346.620 

44414.620 

1171.592 

1183.980 

-12.410 


The computed geoid heights and the positions of the points were applied in equations (2) and (3) to obtain the model parameters 
of Models 1 and 2. While the parameters of model 3, were computed with the geoid heights, positions and the ellipsoidal heights 
of the points using equation (4). The model parameters of the three geometric geoid surfaces were computed using the least 
squares technique, as well as detailed in Eteje and Oduyebo (2018) and Eteje, et al, (2018). Table 2 shows the computed models’ 
parameters. 
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_ Table 2: Computed Models 1, 2 and 3 Parameters _ 

Model 1 Parameters Model 2 Parameters Model 3 Parameters 


ao 

-1.6214917973 

ao 

-72.6223551631 

ao 

-0.8538128360 

ai 

-0.0000250195 

ai 

0.0001301207 

ai 

-0.0000265750 

a2 

0.0000159308 

a2 

0.0018495211 

a2 

0.0000150071 



a3 

-0.0000000040 

a3 

-0.0007055588 


The computed parameters and the coordinates of the points 
were used to develop Microsoft Excel programs using 
equations (2) and (3) (model 1 and 2 respectively). A Microsoft 
Excel program was also developed using the computed 
parameters, coordinates and the ellipsoidal heights of the points 
with equation (4). The Root Mean Square Errors (RMSEs), as 
well as the accuracy of the three geometric geoid models, were 
computed using equation (5). All the computations were done 
using the developed Microsoft Excel programs. The contour 
maps of the three models and the known orthometric heights of 
the points were plotted with Surfer 11 software using the 
Kriging interpolation method. 

RESULTS AND DISCUSSION 

Table 3 presents the three geometric geoid models’ orthometric 


heights and the known orthometric heights of the points. The 
models’ orthometric heights, were obtained by finding the 
differences between the ellipsoidal and the models’ geoid 
height. They were computed to show the minimum and the 
maximum orthometric heights of the three geometric geoid 
models. It can be seen in Table 3 that the minimum and the 
maximum orthometric heights of the three geometric geoid 
models, Model 1, Model 2 and Model 3 are respectively 
1149.954m and 1267.965m, 1149.987m and 1267.987m, and 
1149.965m and 1267.978m. It implies that orthometric heights 
can be respectively obtained using the three models (Models 1, 
2 and 3) within the ranges of 1149.954m to 1267.965m, 
1149.987m to 1267.987m, and 1149.965m t 1267.978m in the 
study area. 


Table 3: Models 1, 2, 3 and Known Orthometric Heights 


STATION 

Model 1 
Orthometric 
Height (m) 

Model 2 
Orthometric 
Height (m) 

Model 3 
Orthometric 
Height (m) 

Known 
Orthometric 
Height (m) 

71Y65 

1222.157 

1222.117 

1222.136 

1222.122 

71Y80 

1265.401 

1265.419 

1265.413 

1265.410 

71Y97 

1267.965 

1267.987 

1267.978 

1267.942 

71Y125 

1149.954 

1149.987 

1149.965 

1150.081 

71Y141 

1188.900 

1188.889 

1188.885 

1188.934 

71Y143 

1165.363 

1165.363 

1165.361 

1165.290 

71Y149 

1164.133 

1164.114 

1164.140 

1164.099 

71Y151 

1249.633 

1249.624 

1249.646 

1249.643 

71Y152 

1223.800 

1223.792 

1223.795 

1223.836 

71Y153 

1205.679 

1205.677 

1205.664 

1205.659 

71Y154 

1186.616 

1186.605 

1186.615 

1186.604 

71Y155 

1183.999 

1184.024 

1184.002 

1183.980 


Figures 3 to 6 present the contour plots of Models 1, 2, 3 and 
the known orthometric heights of the points. They were plotted 
to present graphically, the shapes, as well as the agreements of 
the three geometric geoid models’ orthometric heights with the 
known orthometric heights. It can be seen in Figures 3 to 6 that 


the shapes of the three geometric geoid models of the study 
area are identical with that of the known orthometric heights. It 
shows the agreement of the three geometric geoid models’ 
orthometric heights with the known orthometric heights of the 
study area. 
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Figure 5: Contour Plot of Model 3 Figure 6: Contour Plot of the Known 

Orthometric Heights Orthometric Heights 

Table 4 presents the computed RMSEs of the three geometric geoid models of the study area. The RMSEs of the three geoid 
models were computed with the 12 points used for the determination of the models to compare the accuracy of the three geoid 
models to determine which of them is most suitable for application in the study area. The accuracy of the model is inversely 
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proportional to the computed value of the RMSE. Therefore, the smaller the computed RMSE value, the higher the accuracy. It 
can be seen in Table 4 that the RMSEs of Models 1, 2 and 3 are respective 0.048m, 0.044m and 0.047m. It implies that Model 2 
is most suitable for application in the study area. 


Table 4: Models 1, 2 and 3 RMSEs 


RMSE 

Station 

Dff. b/w Known H 

Dff. b/w Known H 

Dff. b/w Known H 


and Model 1 H 

and Model 2 H 

and Model 3 H 


Squared (m) 

Squared (m) 

Squared (m) 

71Y65 

0.001190754 

0.000020820 

0.000193387 

71Y80 

0.000075383 

0.000085953 

0.000007216 

71Y97 

0.000529574 

0.002002292 

0.001296752 

71Y125 

0.016237504 

0.008786065 

0.013396262 

71Y141 

0.001130201 

0.002005987 

0.002407356 

71Y143 

0.005348817 

0.005367070 

0.005004249 

71Y149 

0.001136665 

0.000232971 

0.001699064 

71Y151 

0.000103819 

0.000377920 

0.000007582 

71Y152 

0.001263377 

0.001916712 

0.001668608 

71Y153 

0.000414214 

0.000332562 

0.000024855 

71Y154 

0.000148775 

0.000002057 

0.000128597 

71Y155 

0.000343770 

0.001944201 

0.000484578 

RMSE (m) = 

0.048237998 

0.043850702 

0.046831708 


Table 5 presents the computed RMSEs of the three geoid models using the three levelling points (stations 71Y121, 71Y126 and 
71Y147). It was also done to determine which of the three models is most suitable in terms of accuracy for application in the 
study area. It can also be seen in Table 5 that the RMSEs of Models 1, 2, and 3 are respectively 0.091m, 0.066m and 0.089m. It 
as well shows that Model 2 is most suitable for application in the study area. 


Table 5: Models 1, 2 and 3 RMSEs Using the Three Points Spirit Levelling Orthometric Heights 


Test Points by Spirit Levelling 

Station 

Dff. b/w Spirit Levelling H 
and Model 1 H Squared 
(m) 

Dff. b/w Spirit Levelling 
H and Model 2 H 
Squared (m) 

Dff. b/w Spirit Levelling 

H and Model 3 H 
Squared (m) 

71Y121 

0.010292824 

0.005761370 

0.012577696 

71Y126 

0.006365151 

0.005570688 

0.004126757 

71Y147 

0.008245191 

0.001609436 

0.007279882 

RMSE (m) = 

0.091110127 

0.065679764 

0.089413526 


Table 6 presents the computed RMSEs of the three geoid models using the four points (U2004 to U2007) whose orthometric 
heights were obtained by gravimetric means. It was as well done to determine which of the three models is most suitable in terms 
of accuracy for application in the study area. It can also be seen in Table 6 that the RMSEs of Models 1, 2, and 3 are respectively 
0.182m, 0.239m and 0.200m which shows that Model 1 with RMSE of 0.182m is most suitable for application in the study area. 


Table 6: Models 1, 2 and 3 RMSEs Using the Four Points Gravimetric Orthometric Heights 


Test Points by Gravimetric Geoid (EIGEN 6C4) 

Station 

Dff. b/w Gravimetric H 
and Model 1 H Squared 
(m) 

Dff. b/w Gravimetric H 
and Model 2 H 
Squared (m) 

Dff. b/w Gravimetric H 
and Model 3 H Squared 
(m) 

U2007 

0.019749066 

0.021835601 

0.022338538 

U2006 

0.036861977 

0.037627664 

0.042998403 

U2004 

0.007899031 

0.158347164 

0.008743013 

U2005 

0.067964632 

0.010797529 

0.085729513 

RMSE (m) = 

0.181985375 

0.239064823 

0.199880882 
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The results, as well as the RMSEs (0.048m, 0.044m and 
0.047m) of the three models presented in Table 4 simply show 
that Model 2 is most suitable among the three models. But as 
the RMSEs of the three models differ with only about 4mm, 
the three models can be applied in the study area with more 
weight attached to model 2. Also, considering that the results 
presented in Table 5 show that model 2 is most suitable for 
application in the study area, the three models can as well be 
applied in the study area for orthometric heights computation 
as their RMSE, as well as accuracy, differ by only about 2.5cm. 
Although, the accuracy of the three models presented in Table 
6 are not as high as those respectively presented in Table 4 and 
Table 5 when the 12 points used for the geoid modelling and 
the 3 levelling points orthometric heights, were used. Also, the 
results presented in Table 6 show that model 1 is the best 
among the three models. It is well known that geoid modelling 
is carried out to replace spirit levelling whose fieldwork is 
tedious and time-consuming (Eteje et al., 2018). So, as the 
results of the spirit levelling validation, as well as the test 
points, agree with those of the 12 points used for the geoid 
modelling, the three models can be applied in the study area 
with more reliability, as well as confidence in model 2. Also, 
the accuracy achieved for the three geoid models agrees with 
those obtained by Oduyebo et al (2019) for Benin City and 
Oluyori et al. (2018) for FCT, Abuja all in Nigeria. The two 
studies applied the same method (geometric method) in small 
areas as well. 

CONCLUSION 

The study has determined the local geometric geoid model of 
Kampala by comparing the accuracy of three plane geometric 
geoid surfaces. The results of the study show that model 2 is 
the best for orthometric heights interpolation in the study area. 
The study has also developed Microsoft Excel programs for the 
application of the models in the study area. The determined 
model is useful to the Ministry of Land, Housing and Urban 
Development (MLHUD), Civil engineers, Surveyors and 
Geophysicists to convert ellipsoidal heights from GNSS 
observation to orthometric heights in the study area. 
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